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Abstract. We examine the equilibrium conditions of a curve in space when a 
local energy penalty is associated with its extrinsic geometrical state characterized 
by its curvature and torsion. To do this we tailor the theory of deformations to 
the Frenet-Serret frame of the curve. The Euler-Lagrange equations describing 
equilibrium are obtained; Noether's theorem is exploited to identify the constants 
of integration of these equations as the Casimirs of the euclidean group in three 
dimensions. While this system appears not to be integrable in general, it is in 
various limits of interest. Let the energy density be given as some function of the 
curvature and torsion, /(K,r). If / is a linear function of either of its arguments 
but otherwise arbitrary, we claim that the first integral associated with rotational 
invariance permits the torsion r to be expressed as the solution of an algebraic 
equation in terms of the bending curvature, k. The first integral associated with 
translational invariance can then be cast as a quadrature for k or for t. 
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1. Introduction 

Consider a curve in space. Suppose that the curve is sufficiently smooth so that 
the Frenet-Serret frame adapted to it is defined. The curvature k and the torsion 
T then provide a complete characterization of the curve; once they are known, it 
can be reconstructed up to euclidean motions. In this paper we examine local 
reparametrization invariant hamiltonians for curves of the form 



where s denotes arclength, and / is any scalar under reparametrizations. 

Such hamiltonians play a role both in the static and in the kinematic description 
of curves. In the former, we interpret the hamiltonian as the energy of the physical 
system; in particular, an energy of the form f = penalizing bending models the 
stiffness of a polymer [Q, || , and it has been used to model the elastic properties of DNA 
(see e.g. ^). When the energy depends only on the curvature, as Bernoulli and 
Euler both knew, the equilibrium conditions are integrable the torsion is always a 
function of k, and k satisfies a quadrature. Our focus will be on the general case (^. 
Whereas a linear dependence on r, associated with a constraint on the total torsion, 
has been considered |^ , very little appears to be known beyond that 0| . At the very 
least, one is interested in a quadratic dependence on both k and r, the second order 
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terms in a Taylor expansion of /. Such terms appear in the hamiltonian describing 
chiral polymers Q. 

An additional motivation for studying these hamiltonians is the role they play 
in the connection between the motion of curves and integrable systems: the equation 
describing the evolution of some function of the curvature and the torsion with respect 
to certain length preserving vector fields coincides with the nonlinear Schrodinger 
equation; other functions give other known integrable equations. These curve 
motions appear in a number of contexts: vortex filaments and patches in fluids 
(H, H^, 0, classical magnetic spin chains |l^, interface dynamics etc. Specific 
hamiltonians of the form ([^) emerge as conserved quantities under these motions 

Our strategy in the study of the hamiltonians will be to exploit Noether's 
theorem to identify the equilibrium conditions as conservation laws associated with 
the euclidean invariance of the energy. These conservation laws in turn permit 
us to identify first integrals of the equilibrium conditions. Remarkably, in certain 
cases, these integrals can be combined to provide a quadrature for either k or r. 
The constants of integration are the Casimir invariants of the euclidean group. We 
show that, in addition to the well known case of a pure bending hamiltonian, a 
pure torsion hamiltonian also leads to integrable equilibrium conditions. This is 
surprising because, in contrast to the curvature which depends on two derivatives 
of the embedding function for the curve, the torsion depends on three. The Euler- 
Lagrange equations which result involve six derivatives; as such one would not expect 
them to be tractable. The torsion is determined by a quadrature. For a polynomial 
/, the potential appearing in this quadrature is a rational function. We identify other 
hamiltonians with a joint dependence on k and t which are reducible to a quadrature. 
In general, unfortunately, it does not appear to be possible to identify a quadrature. 
The integrals of the motion can, however, be used to reduce the equilibrium conditions 
to the motion of a fictitious particle in two dimensions. In any case, this reduction 
should be helpful for studying these systems. 

For most physically realistic materials the local arclength will be constant. 
This is because there will be a large energy penalty associated with stretching the 
curve. Suppose that an arbitrary deformation is decomposed into tangential and 
normal parts. The constraint on the arclength can then be phrased in terms of 
the corresponding response of the tangential deformation to its normal counterpart. 
However, as we will discuss below, a tangential deformation is a reparametrization so 
that the corresponding change in the hamiltonian can always absorbed in a divergence; 
as such, it cannot affect the Euler-Lagrange equations. Thus, whether or not we 
decide to implement a constraint on arclength, the equations themselves describing 
equilibrium are unchanged. 

The paper is organized as follows. In Sect. 2, we begin by giving a self- 
contained account of the theory of deformations of a curve tailored to the Frenet-Serret 
frame. In distinction to earlier work, it is not necessary to implement the constraint 
associated with locally arclength preserving deformations. We obtain directly simple 
expressions for the variation of the curvature and torsion. In Sect. 3, we analyze the 
consequences of the invariance of the hamiltonian under reparametrizations as well as 
under euclidean motions. We show how to obtain expressions for the internal forces 
and torques on any segment of a curve and their relationship with the equilibrium 
conditions. Systems which depend at most on the curvature k is the subject of Sect. 
4. This is extended to systems that depend on the torsion r in Sect. 5. In Sect. 6, 
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systems that depend both on curvature and torsion are considered. In Sect. 7 we 
briefly consider perturbations of the equihbrium conditions. Sect. 8 relates some of 
the results of this paper to recursion schemes that appear in the kinematics of curves. 

2. Curve deformations 

In this section we describe the geometry of embedded curves in three dimensional space 
in terms of the Frenet-Serret basis for the curve, and the effect of a small deformation 
of the curve on its geometry. 

Consider a curve in space described by the embedding x = X(s), where X = 
{X^, X^, X^). The unit tangent to the curve is given by t = X', where prime denotes 
a derivative with respect to arclength s. Clearly, the 'acceleration' t' is orthogonal to 
t. However, t" is not. The classical Frenet-Serret equations 

t' = Kill , 

ni' = - Kt + rn2 , (2) 
na' = - mi , 

describe the construction of an orthonormal basis {t,ni,n2} along the curve. We 
choose an orientation with n2 = t x ni. k and r are respectively the geodesic 
curvature and torsion. The fundamental theorem for curves tells us that the Frenet- 
Serret curvatures k and r determine the curve up to rigid motions [ p7[ . The actual 
curve can always be reconstructed from its curvatures. Thus, they provide a natural 
set of auxiliary variables. Any local geometrical scalar defined along the curve can in 
principle always be expressed as a function of the curvatures and their derivatives. 

We now analyze the change in the geometry of the curve due to an infinitesimal 
deformation of its embedding in space, X(s) X(s) + (5X(s). Let us first decompose 
the deformation into its tangential and normal parts with respect to the basis 
{t,ni,n2}, 

5X = ^fii t -I- ni + 'Sana . (3) 

This is a convenient strategy when one is interested in the variation of 
reparametrization invariant geometrical quantities. Tangential deformations are 
reparametrizations of the curve. We will use the following two facts: the tangential 
deformations of any scalar / and of the infinitesimal arclength are given by 

^ii/^^-ii/', (5||ds = ds^'||'. (4) 

Now consider the normal part of the deformation. A normal deformation of ds is 

S±ds = -ds K . (5) 
This result implies that, for any scalar /, 

5^(/')= + ((5^/)'. (6) 

In particular, for the three scalar functions X, this implies 

5±t = K^-it + (^-irii + «'2n2)' . (7) 

We now use the Frenet-Serret equations to cast as a normal vector (a unit vector 
and its variation are orthogonal): 

d±t^{^i' -T^2)ni + {'^2' + T'^i)n2. (8) 
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Similarly we have 

Si_t' = Kfit' + (^^t)' ; (9) 
using the Frenet-Serret equations and (jsj) we obtain 
6±_t' = - k{^i' -T^2)t 

+ + (k^ - T^)*! - 2t*2' - t'*2] ni 

+ (2r*i' + T'*i + *2"-r2*2)n2. (10) 

We are now in a position to evaluate the normal variations of the two Frenet- 
Serret curvatures. To evaluate S±k, we take a variation of the first of (|^): S±t' — 
{6±K)ni + KS±ni. Dotting with ni we obtain S±k = ni • S±t' . From ( |l0|) we read off 
that 

6^K = + (k^ - T^) - 2r*2' - t'*2 • (11) 

For a planar curve, t = and S±k = 5*1" + k^^Pi. Deformations lifting the curve off 
the plane do not affect the value of k to first order in the deformation. 

To evaluate S±t, we take a variation of the FS equation for n'l and dot with 02. 
We have 

S±T ~ nn2 ■ 5\_t + n2 • 5i_ni . (12) 
We now rewrite the second term on the right hand side as 
n2 ■ 5±_ni = K*i(n2 • n/) + n2 • ((5j_ni)' 
= KT*i + (n2 ■ (5±ni)' 
^1 



-n2 ■ 5i_t' 

K 



(13) 



where we have applied (g|) to n[ , and we have used the FS equations for both n'l and 
n'j. Substituting for n2 • 6±t and n2 ■ S±t', we obtain 

^iT = K(*2' + 2r*i) + |i [2r*i' + r'*i + *2"-T2*2]| • (14) 

For an initially planar curve, d±T = k4'2' + {'^2" /k)' . Only the deformation along the 
direction normal to the plane contributes. Suppose that S±t = 0, then the equation 
K'i> 2' + {'^ 2" / k)' = should not admit any solutions other than those which correspond 
to an euclidean motion of the planar curve. To show this we note that, for a planar 
curve, K = 0', where O is the angle which the tangent makes with, say, the x-axis. 
Then this equation can be recast as 9q4'2 + 'I'2 = 0, with independent solutions, 
^"2 = sin @, cos 9 and \E'2 a constant — a rotation about x, y and a translation. If 
6±T is constant, the only solution is the helix ^'2 oc generated by the planar curve. 
For completeness, let us note that the normals vary according to 

5^111 = - (^-i' - T*2)t + - (1'2" - t2*2 + 2r«'i' + r'^-i) n2 , (15) 

K 

S±n2 = - (^-2' + T^i)t - - (*2" - r^^-a + 2t^i' + t'^i) m . (16) 

K 

So far we have considered arbitrary deformations of the curve. There are special 
deformations that will be of interest in the following. In particular a deformation 
Y = (5X preserves locally the arclength if Jy^s = 0. In terms of the components this 
translates to 

t ■ Y' = - kYi = 0, (17) 
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which imphes that there exists a (non-unique) vector Z such that 

txZ = Y'. (18) 

This is the starting point of the filament model recursion scheme []l6| , where a family 
of locally arclength preserving vector fields {Y(„-)} is defined by t x Y(„) = , 
with Y(o) — — t. We will have something to say about this below in Sect. ||. 

Let us also note that if ((5^X)' • rii = 0, or in components, ^'i' — t5'2 = and 
^2' + T"^i — 0, then S±ni = and d±t = 0; the Frenet-Serret basis is left unchanged 
by this type of deformation. 



3. Invariance and symmetry 

The hamiltonian H for the curve depends locally on the geometry, and it possesses 
various symmetries, both local and global. The local symmetry is reparametrization 
invariance, and it restricts severely the form of H . The global symmetries are euclidean 
motions: translations and rotations. They give rise to conservation laws. 



3.1. Reparametrization invariance 

The hamiltonian H is, in general, a sum of terms, H = Hi + H2 + ■ ■ ■, each of which 
is invariant under reparametrizations of the curve. This results in the form 

H= f ds f{K,T,K',T',---), (19) 



where / is a scalar under reparametrizations, constructed out of the geometrical 
quantities that characterize the curve: the two curvatures, and their derivatives. The 
lowest order in derivatives of the embedding functions non-trivial geometrical models 
depend only on the scalar n. A simple model that penalizes bending of the curve is 
fi=t'-t'~K^. At the next order, a dependence on r as well as k' is admitted. A 
term in / of the form penalizes the torsion of the curve; one of the form (k')^ is 
a higher order differential bending energy and will not be considered. We note that 
the natural hamiltonians /2 = n'l ■ n'l and /a = rij • n'j are given by /2 = + and 
T^. Whereas /i has been considered in considerable detail ||l^, 19 , neither /2 nor 
appear to have been considered. 

The assumption that / is a scalar under reparametrizations implies that the first 
variation of the energy can always be written as 

6H = J ds£,^,+ J ds Q' , (20) 

where Si denotes the normal projection of the Euler-Lagrange derivative of /, and Q is 
the Noether charge {i,j, ■ = 1, 2). The specific form of the first term follows from the 
fact that the tangential variation contributes only in a divergence. Indeed, using (^, 
the tangential part of the variation of the energy is always a total derivative 

Jdsif^^^y. (21) 

This implies that the Noether charge Q, as a linear differential operator which operates 
on the deformation (5X, is of the form 

Q^f^\\ + 2(0)' + */ + •••, (22) 

where, to construct the Q(n)% we use integration by parts to collect in a total derivative 
the normal deformations ^'i and their derivatives. 
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3.2. Translational invariance 

The hamiltonian H is also invariant under the rigid eucHdean motions: translations 
and rotations. Noether's theorem can then be exploited to determine the conditions 
of static equilibrium. Under an infinitesimal constant translation, (5X — e, the energy 
associated with H stored on the segment AB (labeled by its endpoints) changes by 
an amount 

SHab = e ■ / ds [£'n,- F'] . (23) 

J A 

Here we introduce the spatial vector F with 

Q = -e-F, (24) 

and it is constructed by specializing Eq. ( p2| ) to the case of a constant deformation. 

With no external forces acting so that £i = 0, we may identify F as the internal 
force pulling or pushing at a point on the curve segment. (With our conventions, F is 
the force from the part with lower to the one with higher value of s.) In general, this 
force will not be tangential. Translational invariance implies SHab ~ 0. Because the 
endpoints A and B are arbitrary we deduce the local balance of forces 

r n, - F' = (25) 

In equilibrium, this implies the conservation law, F' = 0, i.e. F is a constant vector 
along the curve. We thus associate a spatial vector with each curve, be it open or 
closed. The squared magnitude of this vector, F^, is the first Casimir of the euclidean 
group. The direction along which F points is often indicated by the symmetry of the 
configuration. 

We have come up with three equations of static equilibrium, whereas we only 
possess two independent Euler-Lagrange equations. One of the former must therefore 
be a kinematical statement, or Bianchi identity, associated with the reparametrization 
invariance of H. Let us examine the three independent projections of the equilibrium 
conditions, Eq. (p5|): we decompose F into parts tangential and normal to the curve, 

F = F|| t + Fi ni +F2 n2. (26) 

Then Eq. ( p5| ) is equivalent to the three equations 

i^l'l - K i^i = , 

Fi' + kF^^-tF2^£i, (27) 
F2' + tFi=£2. 

Comparison with ( p7| ) shows that F, seen as a deformation of the curve, preserves 
locally arclength, as expected. The first condition is independent of the Euler- 
Lagrange equations. It is the promised Bianchi identity associated with the 
reparametrization invariance of H. Note that these equations could also have been 
derived directly by considering the balance of forces, as in pof . 

There is a non-trivial intcgrability condition on closed curves associated with the 
conservation law (p5|). Taking its projection onto X in equilibrium, we can immediately 
deduce that 



dsF\\^0. (28) 
on any closed loop. We will comment below in Sect. 3.4. on its geometrical origin. 
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3.3. Rotational invariance 

Under an infinitesimal rotation (5X = 11 x X, we have that the energy H on a segment 
AB of the curve changes by 

SHab = n- ds [£' n,xX- M'] , (29) 

J A 

where the spatial vector M is defined by 

Q = -n-M, (30) 

and it is obtained from Eq. ( p2[ ) by specializing it to the case of a constant rotation. 
We identify M as the torque with respect to the origin acting at a point on the curve 
segment. Rotational invariance of Hab implies 

M' = rn, xX. (31) 

We decompose M into the sum of the couple of F about the origin, plus an intrinsic, 
translationally invariant part, 

M = XxF + T. (32) 

Then the differential torque T satisfies 

T' = F x t . (33) 

We emphasize that this equation does not depend on whether the curve is in 
equilibrium or not. Whereas M is conserved in equilibrium, neither X x F nor T 
is. It is clear, however, that the projection of T onto F, the second Casimir of the 
euclidean group, is conserved 

J = T • F (34) 

where F = F/F, and for future convenience, we choose to fold |F| into the definition. 

As we did earlier for F, we can also decompose T into parts tangential and normal 
to the curve, 

T ^ Til t + Ti ni + T2 n2 . (35) 
Then Eq. ( p3| ) is equivalent to 

T|| - kTi = , 

Ti' - TT2 + kTii = F2 , (36) 
T2' + rTi = -Fi , 

where we are using the convention t • rii x n2 = 1. The first equation, as was the case 
for F in Eq. (p8|), is also valid off-shell, and it says that T, seen as a deformation of 
the curve preserves locally arclength. 

If a curve is deformed along T, i.e. 5'K = T, we find that the corresponding 
variation of the curvature and torsion satisfy 

5Tn = £2, Stt = -(^) . (37) 



Therefore, in equilibrium, not only deformations along F, M corrrespond to rigid 
motions which leave the geometry unchanged, as expected, but also deformations 
along T. 
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3.4- Adapted cylindrical coordinates 



The two conserved vectors F and M together single out a cyhndrical polar coordinate 
system {p,6,z} |2|, |l^, Suppose F ^ 0. We ahgn our coordinate system such 
that F points along the positive z direction. We next perform a translation orthogonal 
to F so that M is rotated into F. Then we have that 

T = JF + FxX. (38) 

Alternatively, we can arrive at this expression by integrating Eq. (p3|), 



j dsT ^ J ds {F xXy - J dsF' xt, 



(39) 



and noting that at equilibrium the second term vanishes, so that T and F x X differ 
by a constant vector. Then contraction with F reproduces Eq. (38), up to a constant 
translation. 

The modulus 



J' 



(40) 



determines the cylindrical radius in the adapted system in terms of and the two 
Casimirs. Typically, will be some function of k and r and their derivatives. 
To complete the construction, we describe the tangent vector in these coordinates, 
t = (p', 0', z') . Thus, the projection t ■ F determines z as a quadrature 

(41) 

This provides the promised geometrical interpretation of the global conservation law 
(p8|). Similarly, the projection t ■ T determines as a quadrature, 

Fp^e' ^ Jz' - T|| . (42) 

The expressions (0), (|l|) and (H) are related by the normalization condition 
p'^ + p^6'^ + z'^ = 1. This is not immediately apparent. It can be shown by taking a 



derivative of Eq. (40) and squaring, together with the squares of Eqs. 
long as p does not vanish. 



4. Bending 

Let us consider an hamiltonian that depends at most on the curvature, 
H = Jds f{K) , 

where / is any local function of its argument. We find that under an arbitrary 
deformation of the curve we have, with = df /dn, 

6H= Jds {fJ^K-fK^ij + Jdsif^^iY 



= d. 



ds 



+ (^2 _^2^f^_^f 



(2t/«)' - 



^'2 

(43) 
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where we have used the expression (|JJ) for S^n. By comparison with Eq. (^0|), we 
immediately read off the Euler-Lagrange equations £i — 0, where 

£i = U" + - r^)f. - (44) 

£2 - i2TfJ - t'/k . (45) 

We see that r contributes to the "driving force" for k in the first equation. Integrating 
the second gives 

/k^t — constant , (46) 

which determines r as a function of k. We show below that the constant appearing 
here is J as defined by Eq. (^4|) . We can substitute into Eq. ( ^ ) for r to obtain a second 
order differential equation for k. It is clear that the Euler-Lagrange equations ( p| ) 
and ( ^5|) are integrable; r is given as a function of k, k is determined as a quadrature. 
The Noether charge Q is identified as the divergence appearing in Eq. (^3|) , 

Q = /*|| + - - 2t/«4'2 . (47) 

The conserved force F is obtained by specializing the deformation to a constant 
infinitesimal translation (5X = e in this expression. In the second term, we use the 
Frenet-Serret equation to obtain 5*1' = e • ni' = e • (— Kt + rn2). Eq. ( ^ ) then gives 

F = - /)t + /k'iii + T/,n2 . (48) 

Note that the tension in the curve is identified as — = / — /kK. It is not constant, 
in general, along the curve, and may also take negative values. 

In a similar way, from the Noether charge corresponding to a rotation (5X = 17 x X 
we obtain the conserved torque, M. In general, only terms with derivatives of the '^i 
contribute to the differential torque T. In this case we have VPi' = (il x X • rii)' = 
fl • (X X rii)' = $7 • (t X rii + X X rii'). The second term contributes to the orbital 
part of M, while from the first we find that the differential torque is given by 

T=-/,txni = -/,n2. (49) 

Note that while the torque due to bending, M, is not generally of the simple form 
F X X unless — 0, i.e. f is constant, neither is T of the most general form: there is 
no component along either t or rii if only bending is penalized. This accords with our 
intuition: the axis of rotation due to the bending which rotates t towards rii is along 
n2. For this model, the second Casimir of the euclidean group, J, given by Eq. (|3^) , 
is read off by dotting Eqs.(|8|) and (H) as 

FJ^-U^T. (50) 

We thus identify the constant appearing in Eq.(^) as —FJ. 

Substituting Eq.(|50|) for r into the magnitude of the force determines 
quadrature, involving the two constants F and J: 

F' = ihr + ih^-fr + ^- (51) 



Suppose /k is not constant. The quadrature (|51|) describes the radial motion (/«) of 
a fictitious particle with a mass — 2, positive energy F^, and angular momentum FJ, 
moving in the central potential, V{k) = {/kK — /)^ . We note that this potential is 
bounded from below. Eq. (|5l] ) can be integrated to determine k implicitly as a function 
of s: 



[F^-v{f,)-F^jyfj]y^ • 



(52) 



Hamiltonians for curves 



10 



Once K, and therefore r via Eq. (]50|), is known, one can use the expressions ([4 Of), 
(pl|), ( ^2[ ) to obtain by a further quadrature the trajectory in the adapted cyhndrical 
coordinates {p,9,z}. In particular, we note that from Eq. ( pO| ) it foUows that the 
radial coordinate p is determined pointwise by the curvature k. 

The physically most relevant model is one quadratic in k, subject to the constraint 
that the total length is constant ||l8| : 

H = J ds{K^+fi), (53) 

where ^ is a Lagrange multiplier enforcing the constraint. We have /(k) = + fi, 
and Eq. (|5^) reduces to 

= 4{.r + - M)^ + ^ , (54) 

where we have used FJ ~ —Ak'^t. The potential becomes quartic for large k. If J 
vanishes, t = 0, and the problem reduces to that of a planar elastica. There is one 
circular solution labelled by the winding number n — ±1, ±2, • • •, correponding to an 
n-fold covering of the circle, with = fi. /i is then fixed by L and n. The vector F 
vanishes on these solutions whereas T does not, pointing out of the plane. The doubly 
degenerate ground states are the once covered circles with n — ±1. 

For n = 0, there are two oppositely oriented figure-eight configurations. There are 
two inflection points where k = 0, connecting equal positive and negative curvature 
lobes. The integrability condition, Eq. (|28|), implies 

ds = fiL , (55) 

where L is the total length of the loop — this fixes fi to be positive. 
Consider next the model described by (see e.g. |^) 

H = J ds{K- Kof , (56) 

where kq is some positive constant, the 'spontaneous' curvature (notice that we do 
not include a constant length constraint in this case). The absolute minimum, H = 0, 
obtains when k = kq, which corresponds to an n-fold circular loop of radius Rq = i^o^ , 
with n arbitrary. The ground state is therefore infinitely degenerate. In this model 
we have, 

F'^4Wr^(>^'-<f^j0^^ (57) 

If J = 0, r vanishes as well, and the potential possesses a single minimum at k = kq. 
The only equilibria with J = as before are the circles and figure eight. If J ^ 0, 
the potential is quite different from the one considered earlier: it now diverges at 
K = kq. On either side of kq, it develops a local miminum. As before, the integrability 
condition can be used to exclude constant k closed loops. Since the 
integrability condition then implies that n — kq if it is constant. 

We end this section with a brief description of the model described by the scale 
invariant bending energy with f = k. In this case, Eq.(^) implies t — 0. We also 
have F = = J. Any plane loop extremizes this bending energy, which (for positive 
k) coincides with the winding number of the loop on this plane. The minimum is 
realized on any convex loop. 
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We turn now to hamiltonians of the form 



H = Jds fir) 



(58) 



where / is an arbitrary local function of its argument. We will see that, like in the 
case H = J dsf{K) discussed in the previous section, such models are integrable by 
quadratures. We will require for the remainder of this section that neither k or r 
vanish. 

We determine the normal variation of the free energy (/r — df /dr): 



6±H\i = 



2-//' + -fr' - 2^fr' -.fn + 2ktU 



*1 



+ / ds 



K K K 



(59) 



SiH\o= ds 



+ / ds 



K 



^^2 



fr 



*2" 



^*2' + (K'-r2):t* 



fr'V 



*5 



(60) 



where we have used Eqs. (||), (0). The Euler-La grange derivatives £i along the normal 
directions are identified as the coefficients of the discarding total derivatives: 



£,=2t{^^] + -fr' - nf + 2KTfr , 



- («/r)' . 



£2 = - { — 



fr'Y , r^fr' 



(61) 



(62) 



As expected, the Euler-Lagrange equations are of order three in derivatives of r, and 
therefore of sixth order in derivatives of the embedding functions. 

The Noether charge is given by the total derivatives in (|59|), (pOh, 



S = /^'ll+2-/,^'i'- 



fr'\' . ffr 



2 - 

K K 



fr 



- hl'2"-— *2' 



This permits us to write down the constant force 



F = [rfr - f)t + -//ni - 



+ K/r 



n2 . 



(63) 



(64) 



We note that the structure of the tangential component is identical to the one of F 
for bending in Eq. (^), with t ^ k. Moreover, the second derivative of r appearing 
in F2 can be lowered to a first derivative by exploiting the Euler-Lagrange equation 
£1 = 0. Note that this was not necessary in first curvature models. We have 
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so that we have the alternative expression: 

F = {rfr - /)t + -//m - ^ 
K 2r 

The differential torque T is given by 



k/ - T 

K 



n2 



-frt 



-ni 



(66) 



(67) 



Notice that in addition to the expected tangential component — /rt due to the twist 
about the 'rod' axis (see e.g. |2^), there is a contribution due to differential twist 
along ni. There is, however, no n2 component. The second Casimir, as defined by 
Eq. (U) , takes the form 



FJ - fr if 



(68) 



Remarkably the solution by quadratures is possible exactly as in the first 
curvature models, as we show now for a special case which is sufficient for our purposes, 



From Eq. 



pw we obtain immediately 



K 

4^2 



Moreover from Eq.(| 
FJ ■ 



K 

we obtain 



(r2 - 2^,f 



(69) 



(70) 



ir(2^ 



T I — 



(71) 
(72) 



which permits to eliminate r'/fc from Eq.(|70[). Doing this, we get 

t'^{t^ - 4/i2) + k^{t^ + Fjf + AFJt'' + 4FV4 = . 
Thus r is determined in terms of J, F and k as a root of an eighth order polynomial. 
It is rather surprising that r is determined pointwise, just as in the pure bending case, 
as some function of k. 

If we insist on adhering to the same mechanical analogue of a non-relativistic 
particle with radial n exploited for bending, the potential appearing in the quadrature 
is going to be a mess. Fortunately, it is also possible to set up a quadrature for t. We 
solve Eq.(p3) for k as a function of r, 



4r^ 



F^ - fi^ + tFJ 



[FJ 



and substitute into Eq.([7l|). We obtain an equation of the form 

r'2 + y(r,i^,J,//)=0, 
where the potential is given by 

(F2_;,2_^rFJ + ir4)(/iT- 



y(T,F, J,m) =4r^ 



FJ 



ir3) 



(73) 
(74) 
(75) 



(FJ + t3)2 

which again describes a non-relativistic particle (this time with position t and zero 
"energy") moving in a potential which is a ratio of polynomials. The analysis of the 
equilibrium configurations for this model is outside the scope of this paper. We note 
that the potential tends asymptotically to r^/2 — a quartic once again. Note that 
the integrability condition 

j ds{T^ - 2^) = (76) 

implies that /z must be positive on a closed loop. 
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6. Bending and torsion 



Let us now consider models with a joint dependence on k and r, / = f{K,T). The 
integrability exhibited by the cases considered so far does not persist, in general, when 
/ depends on both k and r. Unfortunately, this is also the case of interest in biophysics 
where, for example, models of the type / = an^ + I3t'^ are studied. 
The Euler-Lagrange equations take the form 

f 

J J T 



f'r 



2r ( 



t'^ + 2kt/^ (77) 



£2 = = 2(t/«)' -T'f^-['-^\ + T^l^ - {nuy 



(78) 



Note that one has to be careful to include only once the term —nf in adding Eqs. 
(Q) and (|6l| ) to obtain £1. The force and the differential torque are 



UrT + /„K - /)t + -{nf,' + Tfr')ni 
K 



fr' 

frt H ni + f^n2 

K 



K 



n2 : 



(79) 



(80) 



The force is obtained by adding Eqs. ^aj), (p^), taking care to include the term —ft 
only once. The differential torque is given by the sum of Eqs. (p9|), (|67|). Note that 
T has components along all directions. It follows that the second Casimir is 



FJ - -frifrT + f.K - f) - + //t) + 



— ) + Hfr - Tf^ 



(81) 



By comparing this expression for FJ with the corresponding one for the model 
depending only on r, it is clear that the same strategy exploited in Sect. 5 to produce 
a quadrature will not work. There are, however, two interesting mixed cases that are 
tractable by quadratures. 

The first possibility is to consider the bending energy constrained to a fixed length 
and torsion (see for a detailed analysis of this model). Thus, we consider the model 
defined by 

/ = i«2 + ar + /i. (82) 

The total torsion, T = J ds t is dimensionless. The fact that makes this model 
a minimal variation with respect to the pure bending models of Sect. 4 is that it 
does not introduce derivatives of t in the equilibrium conditions. The force and the 
differential torque are given by 



F = i(K^ - 2^)t + k'ui + k{t - a)n2 , 
T = — at — Kn2 . 
We have then that the second Casimir is 

FJ = afi + —{a - 2t) , 
This invariant can be inverted for r as 



1 



r = -^{afi — FJ) 



a 
2 



(83) 
(84) 

(85) 

(86) 
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The corresponding quadrature takes then the form 



.'2 



1 



1 

4k2 



afi — FJ 



(87) 



Note that the torsion constraint does not affect the integrabihty condition Eq. (|28| ' 
ahhough it affects the form of the coordinates p, 9 as defined in Sect. 3.4 
The second possibihty is given by adding a term hnear in k to the model ( |69|) 

that 



so 



■' 2 



fi. 



The two Casimirs now take the form 



F2 = 



FJ 



4r2 
1 



— ~ fi + a — 
2 K 



{2p 



27/ 



an 
27 



T 

K 



(89) 



(90) 



Note that the latter is considerably more complicated than in the pure torsion case, 
as given by Eq. (|7l|). It is possible to use FJ to eliminate t'/k in F^, so that r is 
determined by k pointwise. However, the resulting expression is quite messy. 

Although it is clear that the general case will not be reducible to a quadrature, 
the use of the Casimirs of the euclidean group allows for significant simplifications 
over a direct approach at the level of the equilibrium conditions. We illustrate this 
fact with an example: let us look at the model 

f^li^'+r'). (91) 



This is known as total curvature 



and it is a natural function of curvature and 



torsion, in the sense that ni' • rii' — i?+t'^. It also appears in Ref. |24] as a conserved 



hamiltonian. From Eqs. 

£i = 2r 

£2= - 



^) , we read off the equilibrium conditions 



2^ 



0. 



- + 



(92) 
(93) 



For the force and differential torque, from Eqs. ([79|), (pOj), we obtain 



f4(.^ 



-(rt. 



2k 



ni - 



n2 



ni 



K 



It follows that the Casimir invariants take the form 

2 



^2 ^ 



FJ 



-) -^i- 



(94) 
(95) 



(96) 



Hamiltonians for curves 



15 



We can eliminate the second derivative in F'^ using the definition of J to provide a 
condition of the form 

TiK,K',T,T',F,J) =0, (97) 

which can be considered as the energy condition for the motion of a fictitious particle 
in two-dimensions. 



7. A note on perturbations 

We have, so far, focused on exact methods in deahng with the Euler-Lagrange 
equations for the elastic models we consider. We digress briefly in this section 
to comment on the perturbative analysis of these equations — a more expanded 
treatment of this material will be presented elsewhere Although outside 

the main focus of this article, we find it instructive to include here an example 
of a complementary approach which, through approximations, allows a complete 
treatment, from the energy functional to the actual embedding that minimizes it. 
Apart from the obvious benefit to intuition, we obtain a non-trivial check of many of 
our formulas by computing the force F and the Casimir FJ and verifying that they 
are constant. We do this in the particular case 

f = K^+T'^+fl\ (98) 

i.e. when the total curvature is penalized, with constrained length. 

Using Eqs. ([ztI), ([zsI), the equilibrium conditions can be satisfied, with / as above, 
by constant (non-zero) curvature and torsion, k = kq, r = tq, provided that 

4 + r's = (99) 

(notice that £2 is identically zero for constant curvature and torsion). The resulting 
space curve is a circular helix. Perturbations would give to the axis of this helix a small 
curvature and torsion, while changing, in general, kq and tq as well. We take then k 
and T to be power series in a small parameter e (related to the above "macroscopic" 
curvature and torsion of the axis of the helix) and read off the resulting Euler-Lagrange 
equations order-by-order in e. The zeroth order result is Eq. ( p9| ) above, while to order 
e we get 

n'l + 2tot[' + + Ton = , t[" - Tok[ - t^t[ = , (100) 

where 

K(s) = Ko + eKi(s) + C'(e2), r(s) = Tq + eri(s) ©(e^) (101) 

and we have set kq — 1. The solutions to ( |100| ) involve constant terms, sines and 
cosines, as well as terms proportional to s and s'^. To reduce the number of parameters 
(five initial conditions as well as /i), we choose to eliminate the terms in s and s^, 
resulting in the constraints 

K[iO) = i2^l^ - l)ro-V{(0) , t['{0) = ro(«:i(0) + rori(O)) . (102) 

This amounts to a restriction to periodic solutions, with period equal to that of the 
unperturbed helix. The following abbreviations will be useful in this section 

ai = /i"^(toKi(0) -f a2Ti(0)) , 02 = 2^^ - 1 , as = /i"^(Ki(0) + roTi(O)) . 

We may furthermore set t{(0) = by a suitable shift in s. The solutions then become 

Ki(s) = —aiTo + a-iOLz cos(/is) , 'T\{s) = ol\— asTo cos(/Lts) .(103) 
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We now determine the corresponding embedding, using the Weierstrass representation 
for the curve (we follow the conventions in ) . We first solve the differential equation 

=g(s)$(s), (104) 

where $(s) is the SU{2) matrix describing the rotation of the Frenet-Serret frame and 

Q{s) = -T(s)eo - K(s)e2 , (105) 

with eo = — ei = — §cr^, 62 = ^fcr^ and cr* the Pauli matrices. The embedding 
then is given by 

x{s) + iy{s) = [ 2af3ds' , S{s) = /" {aa - (3i3)ds' , (106) 
Jo Jo 

where a = $11 and (3 = $12. This gives us the helix with its tangent vector, at s = 0, 

along z. To get instead its axis, at s = 0, along z, we rotate around the x-axis by an 

angle rj, with tan 77 — tq/kq. Denoting the resulting embedding by {x, y, z), we find 



x{s) 



fj, ^ cos(^s) + e ^/i ^roQ;i(cos(/is) — 1) — -a3(cos(2/is) — 1) + -TqQ;3S^^ 



(s) = — /i ^ sin(/is) + e ( /i ^Toai sin(/is) — -aa sin(2/is) + (-/ias — /i ^Toai)s 



z{s) — fi Tos + eyn (ai + Tgaa) sin(/is) + /i Toa^ s cos{fis) + fi ai s j . (107) 

We notice that the axis of the helix is bent in the x-z plane due to the term in x{s), 
while the term linear in s in y(s) gives it torsion as well. As a check of our general 
formulas in the previous sections, as well as of our explicit calculations in this one, we 
compute now t, ni and n2 and then the force F from ([79|). We find 

F = ~2eAi3c3y, (108) 

showing that the (indeed constant) force only appears as a result of the perturbation, 
an exclusive feature of this particular model. Finally, Eq. (^ij) gives that FJ is zero 
to this order in e, J = + O(e^). 

A final remark is due concerning the validity of the above first-order solution. 
When moving on to second order, the solutions for K2, T2, will again involve periodic 
as well as non-periodic terms. Eliminating the latter fixes some of the parameters that 
are arbitrary in the first-order results. More generally, requiring periodicity at order 
k, restricts the solutions found to orders less than k, a feature that can be traced to 
the non-linearity of the Euler-Lagrange equations. 



8. Filament model recursion scheme 



In this section, we discuss briefly the filament model recursion scheme, and its 
relationship with the Noether currents for curves. It is described in great detail by 
Langer in the nice review [p^ , and our brief discussion is contained in his work. We 
are only changing the point of view by putting hamiltonians to the forefront, rather 
than curve motions. 

As we mentioned briefly at the end of Sect. 2, consider spatial vector fields Y 
which locally preserve arclength, then the filament model recursion scheme is defined 

by 

t X Y(„) = Y(„_i)' , with Y(o)=-t. (109) 
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The first few terms in tliis liierarciiy are 

Y(2) = yt + K ni + KTn2 , 

^3 

Y(3) = K^Tt + {2tk' + KT')ni + (kt^ - k" - • 

These vector fields have remarkable properties. First, we recognize that Y(i), known 
as the filament model, is the differential torque T for the model / = (1/2)k^, and 
Y(2) is both the force F for the same model and also the differential torque for the 
model / = K^T. Moreover, the integral of the tangential component of the n-th gives 
the corresponding conserved hamiltonian. 

Now, is it possible to set up alternative recursion schemes? From the point 
of view of the hamiltonians, one can start start with some H = [ds fi^K^r), and 
compute its differential torque T set T — Z(i), then from Eq. ( p3|) it follows that 
the associated force is F = '^(2)- Now from the equilibrium conditions in the form 
(p5|), we have f^n^ = t x Z(3), so that Z(3) = ■^(3)||t + £2111 — £in2. However, 
to satisfy the condition that the vector be arclength preserving, we need to satisfy 
Z\\ — k£2 = 0, and this 'perfect derivative phenomenon' appears to happen only in 
the filament model recursion scheme. For example, for the model quadratic in t, k£2 
is not a total derivative. 
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